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ABSTRACT 

Context. Twisted coronal loops should be ubiquitous in the solar corona. Twisted magnetic fields contain excess magnetic 
energy, which can be released during magnetic reconnection, causing solar flares. 

Aims. The aim of this work is to investigate magnetic reconnection, and particle acceleration and transport in kink- 
unstable twisted coronal loops, with a focus on the effects of resistivity, loop geometry and atmospheric stratification. 
Another aim is to perform forward-modelling of bremsstrahlung emission and determine the structure of hard X-ray 
sources. 

Methods. We use a combination of magnetohydrodynamic (MHD) and test-particle methods. First, the evolution of 
the kinking coronal loop is considered using resistive MHD model, incorporating atmospheric stratification and loop 
curvature. Then, the obtained electric and magnetic fields and density distributions are used to calculate electron and 
proton trajectories using a guiding-centre approximation, taking into account Coulomb collisions. 

Results. It is shown that electric fields in twisted coronal loops can effectively accelerate protons and electrons to 
energies up to 10 MeV. High-energy particles have hard, nearly power-law energy spectra. The volume occupied by 
high-energy particles demonstrates radial expansion, which results in the expansion of the visible hard X-ray loop and 
a gradual increase in hard X-ray footpoint area. Synthesised hard X-ray emission reveals strong footpoint sources and 
the extended coronal source, whose intensity strongly depends on the coronal loop density. 
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1. Introduction 

Twisted magnetic fields should be ubiquitous in the so¬ 
lar corona. Various types of photospheric motions, such as 
shear and r otation ( often observed in sunspots, see for ex¬ 
ample Brown et al. 2003) add helicity to the coronal held 
structure. At the same time, the new emerging magnetic 
flux requires a non-zero helicity in order to form a loop, 
ins tead of a diffuse magn etic held arcade (see simulations 
by Archontis et al. 2010 and discussion therein). Twisted 
coronal loops are oft en directly observed in extreme ultra¬ 
violet (EUV) (see e.g. Raouah 2009 iSrivastava et al. 2010) 
and found in reconstructed coronal magnetic fields (e.g. 


Regnier and Amari 2004 Malanushenko et al. 2011). 

The twisted magnetic field contains excess magnetic en¬ 

ergy which can be released with the held relaxing towards a 
potential conhguration. For this reason twisted huxtubes in 
the corona are often associated with solar hares and other 
active events. A number of observational works reveal the 
presence of a twisted magnetic held in pre-hare active re¬ 
gions; some of them directly associate th e twisted magnetic 


ropes with the energy release (see e.g. Srivastava et al. 


2010). Furthermore, there seems to be a strong link between 


Kliem 


dence and simulations (see Amari et al. 2000 Toeroek and 
2005 Gibson and Fan 2006 Fanf 2010| ). However, 


reconnection in twisted magnetic loops m the corona could 
also explain event s with no plasmoid ejection or CME. 
Aschwanden et al. (2009) showed that a substantial num¬ 
ber of small hares occur in a simple single loop conhgura¬ 
tion, also known as self-contained hares. It is suggested that 
magnetic reconnection in twisted coronal loops in an active 
region with no open held could explain these phenomena. 

There are several studies concerning magnetohydrody¬ 
namic (MHD) properties of twisted coronal loops. Kliem, 
Toeroek, and colleagues investigated kink and torus insta¬ 


bilities in twisted loops in zero-beta plasma (Toeroek and 


Kliem 

2005 

Kliem et al. 

2000 

2010 

). The model was 

later applied 

to study the ejective eruptions and CME. 

Browning and Van der Linden 

(2003 

) and 

Browning et al. 


a semi-analytical approach, determining criteria for ideal 
kink instability and subsequent relaxation during a non- 


kink and torus instabilities in the twisted coronal ropes, and 
the ejection of plasmoid and a subsequent coronal mass 
ejection (CME), which is confirmed by observational evi- 


cently generalised by Bareford et al. (2013). These criteria 

has been used by Browning et al. (2008), 
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O 

(2009). 

and later by 

Gordovskyy and Browning 

(20111) to perform 


e-mail: mykola. gordovskyy@manchester .ac.uk 


numerical simulations of magnetic reconnection in a kink 
unstable magnetic huxtube. All these studies consider in¬ 
stability and magnetic reconnection in a huxtube, which is 
cylindrically symmetric. 
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Height, Mm 



Fig. 1. Initial density (panel a), temperature (panel b), and column density (panel c) distributions. Solid lines correspond 
to the low-density case (Model A); dotted lines correspond to the high-density case (Model B). 


Gordovskyy and Browning (2011 2012) show that the 
energy release during reconnection in a twisted loop could 
provide a more convenient model of particle acceleration 
when compared to the standard model. The main benefit 
is that particle energization is distributed within the loop, 
rather than localized in the upper corona and hence, the 
energy losses due to various transport effects (collisions, 
return current) are reduced. Indeed, it is easier to acceler¬ 
ate a larger number of particles in a large-scale distributed 
region rather than in a compact volume, and hence, this 
scenario may have important consequences for the particle 


numbe r problem (see e.g. Brown 1976). Gordovskyy et al. 
(2012) also demonstrated that reconnection in twisted loops 
with magnetic Held convergence near footpoints could pro¬ 
vide a configuration for par ticle re-acceleration in the lower 
corona and chromosphere (Brown et al. 2009). They also 
show that the developed model can explain some features of 
the variation in hard X-ray sources recently observed with 


RHESSI dKontar et al.[ 2011). 


Gordovskyy and Browning (2012) used the model devel- 


ky; _ 

oped in Hood et al. (|2009 ) and Gordovskyy and Browning 


(2011) to study particle acceleration in twisted loops using 
the test-parti cle approach. Th i s stud y recently has been 
extended by Gordovskyy et al. (2012) to incorporate colli¬ 
sions and magnetic field convergence. The obtained time- 
dependent particle distributions are used to predict the 
structure of hard X-ray sources. 

In the present paper we develop a model of energy re¬ 
lease in an unstable twisted coronal loop, focusing both 
on thermal and non-thermal components. The model com¬ 
bines 3D MHD simulations and test-particle calculations 
of electron and proton trajectories in presence of Coulomb 
collisions. The key new feature of this model are the large- 
scale curvature of the twisted magnetic fluxtube, and strat¬ 
ification of the ambient atmosphere. The latter affects the 
MHD part of the simulations through density-dependent 
anomalous resistivity, as well as particle acceleration due 
to Coulomb scattering of high-energy electrons and, to a 
lesser extent, protons. In the present paper we also calcu¬ 
late hard X-ray bremsstrahlung emission, enabling direct 
comparison of our model with observations. 


2. Evolution of kinking twisted coronal loops 

In this section, we consider the evolution of a curved twisted 
coronal loop using a MHD approach. The method and 
model set-up are described in Sect. |2.1[ The derivation of 


the configuration with twisted fluxtube using ideal MHD 
is described in Sect. |2.2[ and the evolution of twisted loop 


after kink instability is described in the Sect. 2.3 


2.1. The model and main equations 

Consider a three-dimensional domain with the dimensions 

— [ ^maxj T^max]: V — [ l/maxi +y max ] and ^ = [0; z max ]. 

The evolution of magnetic field B, electric field E, current 
density j, plasma velocity v, density p and specific thermal 
energy w with time t can be described by the standard set 
of resistive MHD equations: 


II 

<a.| 

cc> 

-V • (pv) 

(i) 

dv 

~dt = 

— (v ■ V)u — —j x B — — Vp + pg 

P P 

(2) 

dw 
~dt = 

-(v ■ V)w - (T - l)w V ■ V + - j 2 + Sother 

p 

(3) 

II 

8I« 

-V X E 

(4) 

E = 

—v x B + ? 7 ji 

(5) 

3 = 

— [V x B} 

(6) 

V = 

To 

(r — l )pw. 

(7) 


Here, /x 0 is the magnetic permeability of the vacuum and 
g is resistivity. The source term S 0 ther accounts for all 
other sources of heating, mostly due to numerical dissi- 
pati on and shock visco sity used in the MHD simulations 
(see Arber et al. 2001). The specific heat ratio is T = 5/3. 
The temperature can be derived from the specific energy as 
T = M(Y — l)/TZ w. Here, TZ is the universal gas constant, 
and M. is the average molar mass equal to M. = 0.5 m p NA, 
corresponding to fully ionised hydrogen. The gravitational 
acceleration, g = — gz , is assumed to be constant in the 
whole domain. The last term in Eq. [3] accounts for all other 
possible sources of heating. In the present simulations, it is 
an artificial viscous dissipation used to stabilise the solu¬ 
tion. 

The set of Eqns m is solved numerically using 
Lagra ngian remap MHD code LARE3D by Ar ber et al.| 
(2001). The numerical grid of 256x256x512 elements 
(along the X-, Y-, and Z-axis, respectively) is uniform in 


2 
































































Gordovskyy et al.: Particle acceleration and transport in twisted loops 


each direction. (The resolution is higher in the vertical di¬ 
rection to reduce the error due to density stratification.) 

The lower boundary of the domain is used as a driver 
to create a conf iguration with a twisted rope, as discussed 
in the Sect. 12.21 Conditions for the four side boundaries are 
based upon the assumption that the plasma is stationary 
(v = 0), and all the remaining parameters are constant 
across a boundary (|f — 0, = 0, 1^=0, where de¬ 

notes a gradient across a boundary). All the parameters (v 
and B , w and p) are constant through the upper boundary. 
In addition, gravitational acceleration is set to z ero at the 
upper boundary to maintain its stability (see e.g. Archontis 


et al. 


2010 ). 

Initially, there is a plain-parallel gravitationally strat¬ 
ified atmosphere. The initial density distribution is con- 
struc ted so that it is simil ar to that in empirical models 
(e.g., Vernazza et al. 1981). At the photospheric level, the 
density is expected to be around 1.67 x 10 -4 kg m~ 3 , and 
the temperature is ~ 6000 K. In the model chromosphere, 
the density decreases exponentially by a factor of 100 for 
every 1 Mm. At z ss 2 — 3 Mm, the temperature drastically 
increases to ~ 1 MK. In the model corona, the scale height 
becomes substantially larger, so that the density decreases 
by a small factor, while the temperature is nearly constant. 
The initial density is defined using the following formula: 


p(z) = pi exp - 


Z ~ Zch 
Z 1 


P 2 exp 


Z - Zch 
Z2 


( 8 ) 


and the corresponding initial pressure can be calculated 
from 


dp 

dz 


= pg- 


(9) 


The temperature is then derived from the gas law. 
Parameters pi, z\, p 2 , and z 2 are chosen so that the ini¬ 
tial density and temperature profiles satisfy the criteria de¬ 
scribed above, while z c h is used to shift the profile in a 
vertical direction. In the present simulations, we use two 
initial atmospheres: low-density with the coronal density 
corresponding to quiet corona (~ 10 9 cm -3 ) (Model A 
thereafter), and high-density with the coronal density cor¬ 
responding to cases of extremely high density above ac- 
tive regions (~ 10 ln cm -3 ) (Model B thereafter) (see e.g. 
Tripathi et al. 2008). The parameters for Eq. |8|are given 
in Table |2.lj and "the density and temperature for Models 
A and B are shown in Fig. [T] In addition, Fig. [l] shows the 
column density defined as 


«■»-/ 


p{z 


-dz. 


( 10 ) 


where m p is the proton mass. 

By introducing the scale length L 0 , characteristic mag¬ 
netic field B 0 and characteristic plasma density po one 
can derive the dimensionless variables for all used physi¬ 
cal parameters. Thus, B* = B B// 1 , j* = j poLqBq 1 


E* = E p 0 
t 


1 / 2 j/ 2 b ~ 2 


Po 


-1/2 -1/2 r _i 

-On 


P = P PO, 


1/2 1 / 2 R _1 
Po Po B 0 


— t BqPq ' Po ^o 


w 

>2 


= w popB 0 2 , and p* = 


PPoBq 2 , and g* = gpoPoLoB q, where starred variables de¬ 
note dimensionless physical parameters. The position vec¬ 
tor and all the length variables are scaled as s* = sL jj -1 . 


The dimensionless resistivity can be derived as p* = 
p p 0 1/2 p 1 0 /2 Lo 1 Bo 1 . Later in the paper we refer to the 
dimensionless parameters (dropping the star symbol) ev¬ 
erywhere, except for the cases where units are given ex¬ 
plicitly. In the present simulations, we use the following 
scaling parameters: Lq = 10 6 m, Bq = 4 x 10 -3 T, and 
Po = 3.3 x 10 -12 kg m -3 , yielding the characteristic veloc¬ 
ity i>o = 1.95 x 10 6 m s _1 and time to = 0.53 s. 

2.2. Derivation of configuration with twisted fluxtube 

A magnetic field configuration with a twisted, bended flux- 
tube is constructed by adding a twist to an initially poten¬ 
tial magnetic field. 

The initial configuration is the potential field of two 
magnetic monopoles located outside the numerical domain: 


B(t = 0) = .BiLg 


r — mi 
|r — mi | 3 


r - m 2 \ 

I r - m 2 \ 3 ) ’ 


( 11 ) 


where r is the position vector, and mi and m 2 are the 
positions of the monopoles: mi = [0; a; —h] and m 2 = 
[0;— a;— h\. The depth of the “monopoles”, h, determines 
the convergence of the magnetic field towards footpoints. 
In our simulations, a = 6.4Lo and h = 0.2 Lq, providing the 
convergence ratio (the magnitude of magnetic field at the 
footpoint to that at the loop-top) of ~ 10. 

The helicity is injected into the system by applying ro¬ 
tation at the footpoints. Then we follow the evolution of 
the magnetic field using ideal MHD. The centres of rota¬ 
tion are the points with the maximum magnetic field at the 
lower boundary: x = 0 and y ss ±6.4Lo- Rotational velocity 
corresponds to the solid rotation case (see Fig. [2]), where 
Vrot ~ r inside the footpoint radius R = 0.5Lq : 


^roi(^) — OJtwistT X 

1 — tanh 


r — R 
0.05 L 0 J 


tt 


■ exp - 


tt 


(12) 


where untwist and tt W i S t are the characteristic angular veloc¬ 
ity and time, respectively. 

The rotation is applied through v x and v y components 
at the lower boundary, while the vertical velocity is zero. 
The magnetic field components are kept constant across the 
lower boundary {dB/dz = 0), along with the specific inter¬ 
nal energy {dw/dz = 0). The density at the lower boundary 
is kept constant, p(z = 0,t) = p(z = 0,f = 0), and similar 
to the upper boundary, the gravitational acceleration is set 
to zero here to maintain boundary stability. In all the exper¬ 
iments, we have used urtwist = 0-1 rad/s and t tW i S t = 250to- 
Hence, the rotational velocity is low, iHwist ~ 10~ 2 ao, and 
the field structure, effectively, goes through a sequence of 
quasi-equilibrium states. After ttwist ~ 10 3 to; the total 
twist is about 87r, and the magnetic configuration is nearly 
force-free but apparently kink-u nstable. The obtained fiel d 
geometry is similar to that in Gordovskyy et al. (12012), 
except we now have a more realistic bended loop ana grav- 
itationally stratified atmosphere. 

The evolution of the models during and after the onset 
of kink instability is considered using resistive MHD, as 
described in the Sect.2.3. 
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Table 1. Parameters for the initial density distribution defined by Eq. [i] 



pi, kg m d 

zi, m 

p 2 , kg m 4 

Z2, m 

z c h, m 

Model A 
Model B 

1.7 x 10" 4 

2.2 x 10 5 

1.0 x 10 ii 

5.0 x 10 -10 

2.0 x 10 6 

1.5 x 10 s 



r / L 0 

Fig. 2. Maximum rotational velocity (solid line) and final 
twist angle (dashed line) as functions of the distance from 
the centre of rotation. 



X 


Fig. 3. Magnetic field lines during kink instability develop¬ 
ment: view from the top at t w 800fo- The blue colour de¬ 
notes the field lines connecting the footpoint centres; green 
colour denotes an ambient field. 


2.3. Magnetic reconnection in twisted coronal loops 

The kink instability occurs approximately at t = 920to and 
t = 950to in the Model A and Model B, respectively. In 
both cases, the m = 1 kink mode dominates, resulting in 
the sigmoid structure (Fig. [3]). 

The simulations are switched from an ideal to resistive 
MHD regime at t = 896fo- The resistivity consists of two 
components: uniform and constant background resistivity 
dback and non-uniform current-driven anomalous resistivity 
r/ cr . The need to specify some functional form for anoma¬ 
lous electric resistivity (or magnetic diffusion coefficient) is 


an intrinsic flaw in MHD simulations of magnetic reconnec¬ 
tion; the resistivity in most numerical models is arbitrary to 
some extent. Many studies, however, attempt to use some 
physically motivated form of anomalous resistivity, most 
often assuming that it i s caused by the development of ion- 
acoustic tur bulence (see Ugai 1992 Uzdensky 2003 Barta 
et al. 2011). In the present simulations, we adopt similar 
approach and use current-driven anomalous resistivity. The 
total resistivity is defined as 


V = Vback + VcrSign(j - j cr ). 


(13) 


The critical current j cr can be derived from the criteria 
vdrift > u er , whi ch is normally used in such simulations 
(e.g. [Kliem et al.|[2000). Here, v r i„ft = j/(en) is the elec¬ 
tron drift velocity and the critical velocity equal to the 
sound speed v cr = 
current can be written as 


r^. Hence, the expression for a critical 



(14) 


and by using the scaling parameters introduced above, this 
formula can be reduced to a more convenient form 


jcr _ R_n To P 

jo V 2 P9 R L ,pPq' 


(15) 


where /3 S = 2 ^°^ T is the gas to magnetic pressure ratio and 

Rl p = the Larmor radius of a thermal proton. 

The problem, however, is that in global MHD models the 
current density is limited by the numerical grid resolution. 
In our simulations the grid resolution is SL = Lq/256 « 
4 x 10~ 3 To or ~ 4 km in real terms, while it is widely 
accepted that a realistic thickness of a RCS is about the 
Larmor radius of the proton R ^ pi which is ~ 1 m in the 
corona. Hence, the grid resolution is by a factor ~ 10 3 —10 4 
larger than realistic RCS thickness, which means that the 
real current densities should be ~ 10 3 — 10 4 higher. To 
overcome this problem, we introduce a correcting factor of 
8L/Rl p = 4x 10 3 for j cr resulting in the following form: 


jcr _ /l( „ To _P_ 

jo “ \2 P °6Lp 0 ' 


(16) 


which brings the critical current just above the current den¬ 
sities before the kink instability. 

The background resistivity is pback = 10~ 6 . This value 
is substantially higher than real classical resistivity in the 
corona (which would be about 10 , if we consider our 

scaling parameters), but it is not expected to affect the 
model, as the corresponding diffusion time is tdi.ff ~ 
(SL) 2 /pback ~~ 10 2 — 10 3 to, where SL is the grid step. As 
far as the value for the anomalous resistivity is concerned, 
it is more difficult to determine a physically-justified r] cr . 
Therefore, the experiment with low-density atmosphere 
(Model A) has been repeated with several values of 771 : 
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Time Time 

Fig. 4. Total magnetic energy and total kinetic energies as functions of time for the low-density (solid lines) and the 
high-density models (dashed lines). Time is measured in to units, while the energies are measured in \no -BqLq. 


10 -4 , 4 x 10 -4 , 10 -3 and 4 x 10 -3 . As in earlier MHD 
experiments with Lare code, we did not find any notice¬ 
able difference in the equilibrium formed after the recon¬ 


nection (Gordovskyy et al. 2011). The main difference is 


the speed of relaxation: it is nearly proportional to rye/ 5 - 
Additionally, it was found that dissipation due to numeri¬ 
cal viscosity increases for higher ?y CT . values. This can be ex¬ 
plained by faster plasma flows in models with faster energy 
release. This effect was found to be especially significant in 
the experiment with r] cr = 4 x 10” 3 : the energy lost due 
to (artificial) viscous dissipation was comparable to Ohmic 
heating. Therefore, for particle trajectory calculations in 
Section 3 we use experiments with rj cr = 10~ 3 . 

Once the held goes unstable, the current density dras¬ 
tically increases near loop-top and footpoints. The current 
density is high in these two regions due to different rea¬ 
sons: near footpoints strong currents appear due to mag¬ 
netic held convergence, while the current density near the 
loop-top becomes high because of the kink resulting in for¬ 
mation of magnetic rope. These regions, however, are not 
compact: they are extended along the loop legs and from 
time to time the footpoint regions get merged with the loop- 
top one. The increase in the current density above critical 
threshold switches on the anomalous resistivity, which, in 
turn, results in magnetic reconnection and dissipation of 
magnetic energy. This is qualitativel y similar to the recon¬ 
nection scenario in previous studies dBrowning et al.[ [2008 


Hood et al. 2009[ |Gordovskyy and Browning, |2012p. 


ments 

(e.g., Hood et al. 

2009 

Gordovskyy and Browning 

2011 

Bareford et al. 

2013), 

the current gradually be- 


t ss 1900fo, apart from the footpoints, where it remains 
below j cr level. 

Total magnetic and kinetic energies for the MHD mod¬ 
els are plotted versus time in Fig. [4] The magnetic energy 
variations are similar in both models. The difference in the 
maximum energy and the energy of the relaxed state is 
due to the density-dependent anomalous resistivity current 


threshold. The latter is higher in the high-density case, 
and this is why the current density and magnetic energy 
can reach a higher value before the reconnection starts. 
Similarly, as the reconnection proceeds, the maximum cur¬ 
rent density in this case gets below the critical current 
threshold earlier, while the magnetic energy is still rela¬ 
tively high. As far as the kinetic energy is concerned, the 
low- and high-density cases are different during the helicity 
injection phase. Both model loops demonstrate weak kink 
oscillations; however the frequency of oscillations and their 
amplitude are substantially smaller in the high-density case 
(Model B) than those in the low-density case (Model A). 
After the kink occurs (at approximately t > 900 fo), the ki¬ 
netic energy variations are similar: they peak at ~ 0.002 of 
initial model energy and then gradually decay. The kinetic 
energy at this stage is associated largely with reconnection 
outflows. 

In contrast with previous simulatio ns of kink instabil¬ 


ity in coronal loops (Kliem et al. 2010), we do not observe 


Just after the kink, the current is concentrated in a 
cylindrical structure following the loop legs. At the top 
of the loop, there is a current ring corresponding to the 
magnetic plasmoid. The current density is particularly high 
near the footpoints (due to magnetic held convergence) and 


separation and emergence of the plasmoid. The difference 
apparently appears because magnetic reconnection in our 
simulations occurs in a completely closed magnetic config¬ 
uration, and the overlying ambient magnetic held lines are 
tied either to lower (photospheric) or side boundaries. The 
reconnection happens mostly at the very top of the kink¬ 
ing loop near the upper part of the plasmoid ring between 
the twisted held lines and the overlying ambient held lines. 
This results in the reduction of a twist and radial expansion 
of the loop, which is similar to the previous simulations of 


the initially cylindrical loop (Gordovskyy and Browning 


2011 2012). The main two differences with these previous 


simulations are: (a) the process now, obviously, is not cylin- 
drically symmetric and the current density varies strongly 
both along the loop and within its cross-section; (b) the 
absence of magnetic diffusion very close to the footpoints 
mean the two footpoints remain connected in the final equi¬ 
librium, and all the held lines remain connected to both 
footpoints (unlike in Gordovskyy and Browning (2011)). 
The latter is explained by the sketch in Fig. [6] 


3. Electrons and protons in a twisted coronal loop 

In this section, we consider particle motion in reconnect¬ 
ing a twisted coronal loop based on the electric and mag- 
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Fig. 5. Selected magnetic field lines (left panels) and current density iso-surfaces (j = 0.5) (right panels) during magnetic 
reconnection in the Model A (low-density case). Different colours are used for magnetic field lines to demonstrate the 
change of connectivity. Blue lines originate almost at the centre of the footpoint at y ss 6.4; red lines originate almost at 
the centre of the footpoint (y —6.4); other lines belonging to the twisted fluxtube are shown in green. The corresponding 
times are shown at the lower left corners. 
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Fig. 6 . Schematic sketch showing change of connectivity during magnetic reconnection in 


the twisted loop. 


netic fields, and pl asma density obtained in the MHD sim¬ 
ulations (see Sect. 2.3). The methodology is explained in 
Sect. 3.1 w hich followed by a descripti on o f the trajec¬ 
tories (Sect. 3.2), ene rgy spectra (Sect. 3.41, and spatial 
distribution (Sect. 3.5). 


3.1. Test-particle motion in presence of Coulomb collisions 

Let us consider the motion of charged test-particles using 
the guiding-centre approximation with the particle trajec¬ 
tory described by the velocity parallel to magnetic field 
U||(f), drift velocity perpendicular to the magnetic field 
u(t), and Larmor gyration velocity v g {t). This approach 
is valid in this model, since the magnetic field is non-zero 
throughout the domain, and the proton and electron gyro- 
radii are much smaller than the characteristic scale Lq. 

The trajectories of changed particles in the presence of 
collisions can be described by the following set of equations: 


dr 

dt 


d im) 

dt 


u + T^b 


u = u E + 


7 

m ( 7 ^|| ) 2 

q ”/k 2 B 


[6 x (b ■ V) 6 ] 


TO p 

q 7 k 2 B 


[bxV(nB)] 


-E-b — — (6 • V(kB)) 
m 7 

6a 


St 



Sv 

~\r 

a — 

coll 

St 


- coll 


dp 

v(l — a 2 ) Sv 


av 2 5a 

dt 

cq 

coll 

B St 


(17) 

(18) 

(19) 

( 20 ) 


- coll 


Here, r is particle position, p is specific magnetic moment 
p = Vg /B , ue is the ExB-drift velocity ue = E x B/B 2 , 
and b is the magnetic field direction b = B/B. Apart from 
the last two terms in Eq. [l9|and Eq.[20j th is is the st andard 
set of relativistic guding-centre equations (Northrop 1963). 
The relativistic factors are defined as 


7 = 


c 2 — vf, — u 2 — 2 pB 


and 


k= Jl-u 2 E /c 2 . 


The collisio nal terms in the equations above are defined 
below by Eqs. (22) and (25). Let us discuss them and their 
numerical implementation in more details. In collisional 
terms, we ignore the relativistic factor 7 , since particles 
with a value of 7 noticeably greater than 1 are not signif¬ 
icantly affected by Coulomb collisions. The average colli¬ 
sional energy loss of a particle with mass to moving with 

speed v > v th = yj 


lk B T 


in fully ionised thermal plasma 


with temperature T can be expressed as (see e.g. Emslie 
1978): 


d£ „ / . to n 

_2ne 4 A— -, 
al m e c 


( 21 ) 


where A is the Coulomb logarythm. For non-relativistic par¬ 
ticles, the kinetic energy £ = mv 2 /2 and Eq. 
rewritten as 


21 


can be 


dv „ m e n 

= -£- 2 

at m 


( 22 ) 


where K. = 4:Tre 4 A/m 2 and A « 20. The E 
ations in the full velocity Sv/St for Eqs. 
full velocity in the right-hand side of Eq. 


J 22 | yields vari- 
"Tand [20l The 


221 is calculated 


as v = Jv m + 2 pB. Here the drift velocity is disregarded, 


since (a) it is normally comparable to or lower than the 
thermal velocity vthi and (b) the main component of the 
drift, ue, is caused predominantly by bulk plasma motion. 

The pitch-angle distribution of high-energy particles in 
thermal plasma changes with time due to Coulomb colli¬ 
sions as 


df m n d 

dt m e v 3 da 



(23) 


where the pitch-angle is defined as a = v\\/v. In terms 
of individual particles, the pitch-angle diffusion represents 
random deflections with respect to a chosen direction. Let 
us evaluate the probability of a test-particle to be deflected. 
Equation [23] means that the fraction of particles changing 
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Fig. 7. Typical electron trajectories, their energies and 
pitch-angles. Particles are injected with thermal energies 
at the locations denoted by squares on the lower panel. 
Blue, green, and red lines denote particles, which are non- 
thermal at some point, while grey lines denote particles, 
which remain thermal throughout simulations. See text for 
Srore details. 







Y, L 

o 

Fig. 8. Typical proton trajectories, their energies and 
pitch-angles. Particles are injected with thermal energies 
at the locations denoted by squares on the lower panel. 
Blue, green, and red lines denote particles, which are non- 
thermal at some point, while grey lines denote particles, 
which remain thermal throughout simulations. See text for 
more details. 
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its pitch-angle from a to a + Sa within the time interval St 
is 


110, a, Aa) = StJC 


m n 1 — (a + Aa/2) 2 
m e v 3 Aa 2 


(24) 


This, in turn, means the following: if individual particles 
are deflected by Aa with probability n(u, a, Aa) as per 
Eq(24j then the pitch-angle distribution of the whole parti¬ 
cle population should satisfy Eqj23| 

In the numerical scheme used here, the pitch-angle dif¬ 
fusion is implemented through stochastic jumps of a par¬ 
ticle pitch-angle by a fixed value of either Aa = 0.05 
or Aa = —0.05 after every timestep St with the proba¬ 
bility of n(u, a, Aa) or n(u,a, —Aa), respectively (which 
means the pitch-angle remains the same with the probabil¬ 
ity 1 — n(t), a, Aa) — n(u, a, — Aa)). Hence, the variation 
in pitch-angle with time can be written as: 


Sa 

— = Aan(i), a, Aa) — Aan(r, a, — Aa) 

St 


(25) 


with probabilities n as per Eq. [24] 

In each numerical experiment, we calculate trajectories 
for ~ 10 6 test electrons and protons. Initially, particles are 
uniformly distributed within the simulation domain, have 
Maxwellian velocity distributions corresponding to the tem¬ 
perature of 1 MK, and are uniformly distributed with re¬ 
spect of the pitch-angle a from —1 to 1. 

The calculations have been performed using the GCA 
code based on the second order Runge-Kutta scheme 


(Gordovskyy et al. 2010, 2011). In addition to the usual 


limitations on the integration timestep (St <C Sr/v and 
St -C v/a, where Sr is the grid step and a is the accelera¬ 
tion), the collisional terms add another requirement: St <C 

, which is necessary for n(u, a, Sa) < 1. 
(When n(u,a, Sa) ~ 1, the scheme used for pitch-angle 
scattering calculations becomes unstable, and, obviously, 
the probability cannot be n(u, a, Sa) > 1.) 

The time-dependent electric and magnetic fields and 
their spatial derivatives for the right-hand side of Eqs. m 
l20l are taken from the resistive MHD simulations described 
in Sect. |2.3| The data for each particle position is derived 
by linear interpolation of the data within four-dimensional 


(a:, y , z, t) cells from the ajacent grid points (see 

Gordovskyy 

et al., 2011: Gordovskyy and Browning, 2011 

Gordovskyy 

et al., 2012). The domain boundaries are closec 

for thermal 


£ < 1 keV and open for higher energy 
particles. Each particle is followed until the end of simula¬ 
tions at t = 1800to or until it leaves the domain through 
one of the boundaries. 


3.2. Particle trajectories 

Typical electron and proton trajectories are shown in the 
Figs. [7] and [8j respectively. In general, both species are very 
adiabatic and hence, behave very similarly, The main dif¬ 
ference, of course, is that protons are much (by factor of 
~ 43) slower. As a result, their trajectories are smoother, 
since the effect of small-scale fluctuations (in the parallel 
electric field value, magnetic field curvature etc) is negligi¬ 
ble due to time averaging. 

The majority of particles remain in the thermal distri¬ 
bution; only a few particles (~ 5%) are accelerated beyond 


1 keV. This is comparable to typical acceleration efficiencies 
derived from observed hard X-ray emission. This validates 
the use of a test-particle approach, as the effect of high- 
energy particles on the magnetic and electric field should 
be considerably low. 

Particles move predominantly along the magnetic field 
lines. Despite the connectivity changes, all the field lines 
of the twisted loop remain connected to both footpoints, 
and there is no open field. As a result, during the reconnec¬ 
tion most of high-energy particles remain in or around the 
twisted loop and precipitate towards one of the footpoints. 


Particles accelerated by a parallel electric field in an 
almost collisionless corona have a very narrow pitch-angle 
distribution around a = ±1, which means they are strongly 


collimated along the magnetic field ( 

Gordovskyy et al. 

2011 

Gordovskyy and Browning 

2012 

. However, they get 
chromosphere, get- 

scattered due to collisions in the 

ienser 


ting a wider pitch-angle distribution. As a result, there is 
a small but noticeable fraction of high-energy particles re¬ 
flected by the converging field back to the corona. Electrons 
with energies up to ~ 100 keV and nearly all the pro¬ 
tons are thermalised before reaching the lower boundary. 
Hence, only a small fraction of energetic electrons (with 
£ > 100 keV) and some particles accelerated near the foot¬ 
points (which did not have enough time to thermalise) can 
get to the photosphere (i.e., the lower boundary of the sim¬ 
ulation domain). 


3.3. Particle energy spectra 

There are two main factors affecting particle energies: par¬ 
allel electric field and Coulomb scattering. The appearance 
of a strong electric field is a transient and local effect, 
while collisional deceleration of high-energy particles is al- 
ways present. Henc e, unlike in t he collisionless model s by 
Gordovskyy et al. (2011) and Gordovskyy and Browning 
( TO T it is not possible to characterise the acceleration 
process by the final spectra in this type of model, as these 
spectra are inevitably thermal. We consider the energy 
spectra at three different stages of the reconnection process: 
just after the kink instability occurs, during the fastest en¬ 
ergy release (i.e., when d£ m /dt is highest), and during the 
decay stage. The energy spectra for protons and electrons 
for low- and high-density cases are shown in Figs. HUMS 
andH2l 

The electron energy distributions at the beginning of 
reconnection are similar to those obtained in simulations 


with no collisions (Gordovskyy and Browning 2011): the 


spectra are combinations of a Maxwellian thermal distri- 
bution and nearly power-law high-energy tail. However, in 
the low-density case the high-energy tail, surprisingly, ap¬ 
pears to be softer than in previous studies; its power-law 
index is about 2.0 — 3.5, which is observed in many flares. 
In the high-density case, the “high-energy” tail is harder, 
the spectral index is about 1.5 — 2.0. At the later stages 
(t = 1100 — 1500 to), the electric fields gradually de¬ 
cay and the collisions become dominant. This results in the 
hardening of the spectra around a few keV, and at some 
point, a gap appears between the thermal part and high- 
energy part. This spectral hardening at lower energies (of 
few deka-keV) is similar to that, which appears in thick- 
target models. 

Comparing electron and proton energy spectra demon¬ 
strates the contrast in acceleration times: electron non- 
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t=928 t=960 






Fig. 9. Electron energy spectra in the Model A. Times are shown above corresponding panels. 
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Fig. 10. Electron energy spectra in the Model B. Times are shown above corresponding panels. 


thermal spectra are formed within ~ 1 — 10 to after a 
kink, while it takes around 10 — 100 to to form a smooth 
high-energy tail for protons. At the same time, protons lose 
their energies due to collisions faster than electrons with 
the same energy (since d£/dt ~ \Jm/£). As a result, it is 
more difficult in a collisional plasma to accelerate protons 
than electrons. The difference in proton and electron ac¬ 
celeration efficiencies can also be explained in terms of the 


electric held required to accelerate particles with energy £ 
in presence of collisions in plasma with density n. Indeed, 
particles are accelerated if eE > 27re 4 A —n/£, or if the 
minimum held required to accelerate particles is 


E C rit > 27re 3 A — n/8. 

m e 
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Fig. 11. Proton energy spectra in the Model A. Times are shown above corresponding panels. 
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Fig. 12. Proton energy spectra in the Model B. Times are shown above corresponding panels. 


(For thermal electrons {£ = ksT, m = m e ), this becomes 
the standard expresion for Dreicer field.) Obviously, this 
Ecrit value for protons is higher by a factor of m p /m e 
than E cr it for electrons. Therefore, when the typical electric 
field in the system is higher than both E cr i t (protons) and 
Ecrit (electrons), the proton acceleration efficiency should 
be similar to that for electrons. However, when the typical 
electric field is lower than E cr i t (protons) but higher than 


Ecrit (electrons), the number of high-energy protons should 
be lower than that of electrons, which is the case in the 
Model B. 

This difference between the low-density and high- 
density cases can also be seen in total energetics. Thus, 
the total energy carried by high-energy protons and elec¬ 
trons in the Model A (low-density case) is nearly equal and 
is about 6 — 8% of the energy released during the recon- 
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nection, (although the maximum total non-thermal proton 
energy is reached slightly later then maximum total non- 
thermal electron energy). At the same time, the total en¬ 
ergy carried by non-thermal electrons in the high-density 
case (Model B) reaches ~ 4% of the energy released during 
reconnection, while protons carry less than 1%. These val¬ 
ues, however, show the energies contained in non-thermal 
particles at a given moment and do not necessarily reflect 
the total energies deposited by non-thermal particles (or, 
in other words, the electromagnetic energy converted into 
thermal energy through non-thermal particles). This would 
depend on the value of the electric field compared to the lo¬ 
cal E crit value. Thus, particles accelerated by electric field 
E ~ E cr n are almost instantly thermalised, which means 
that the number of non-thermal particles in the system at 
any time remain low, although a substantial amount of en¬ 
ergy is converted from electric to thermal in this process. 

Here, we define non-thermal as particles with energies 
> 1 keV (the temperature of the initial distribution corre¬ 
sponds to around 100 eV). The energy spectrum below this 


energy is not correct in the present model, since Eq. 21 


is 


only valid for non-thermal electrons with an energy that is 
much higher than ksT for the ambient plasma. 


3.4. Particle pitch-angle distributions 

In collisionless models of particle acceleration by paral¬ 
lel electric fields, the pitch-angle distribution is trivial: 
since perpendicular velocities remain ~ vth and paral¬ 
lel velocitites increase drastically, all non-thermal parti¬ 
cles are strongly collimated along magnetic field lines (i.e., 
Vg/v || << 1), apart from particles being reflected by con¬ 
verging magnetic field. Collisions, however, result in effec¬ 
tive scattering of supra-thermal particles with moderate en¬ 
ergies at ~ 1 — 10 keV. Pitch-angle distributions are shown 
in Figure [13] 

I t c an be seen that particles in the low-density case 
(FigH a-b) are also well-collimated along the field lines, 
although the collimation is not as strong as in collision¬ 
less models. The majority of particles (~ 80%) have pitch- 
angles |a| > 0.75. In the high-density case, however, pitch- 
angle distributions are more isotropic: ~ 40 — 50% of par¬ 
ticles have pitch-angles of |a| < 0.75. 

The pitch-angle distributions also reveal the preferred 
directions of acceleration for diferent species. Thus, protons 
tend to have positive pitch-angles indicating that they move 
predominantly towards y = +6.4Lo footpoint, while elec¬ 
trons tend to have negative pitch-angles and, hence, move 
towards the opposite footpoint. This effect can also be seen 
in spatial distribution (see next section). Coulomb collisions 
strongly affect particles in Model B, and, as the result, the 
proton-electron footpoint asymmetry almost disappears in 
the high-density case. 

3.5. Particle spatial distribution and synthetic hard X-ray 
intensities 

The helicity does not change its sign in the initial configu¬ 
ration, which means the current has preferential direction. 
This results in charge separation at the early stages of re¬ 
connection: electrons and protons are accelerated towards 
different footpoints, producing the electron-proton asym¬ 
metry of footpoints that can be seen in Figure [14] As the 


Electrons 



-10 -50 5 10 

Y 

Protons 



Fig. 14. High-energy {£ > 5 keV) electron and proton dis¬ 
tribution in X — Y plane. Most of the particles are concen¬ 
trated at low heights close to the footpoints. 


current structure becomes filamentary with time, the accel¬ 
eration process becomes more chaotic and footpoint asym¬ 
metry decreases. At the end of reconnection, only a slight 
assymetry can be seen. This is consistent with observations 
showing time delay in appearance of a second footpoint 
source. 

Figure [l4] also shows that the footpoint area (which rep¬ 
resents the cross-section of the volume occupied by high- 
energy particles just above the photosphere) increases with 
time. During the reconnection, the radius of footpoints in¬ 
creases from ~ 0.5 Mm to nearly 1.2 Mm. This effect is 
mainly due to the reconnection between the field lines of 
the twisted fluxtube and ambi ent fie ld lines (see discussion 


in Gordovskyy and Browning 2011). 

Based on the electron distributions derived in Models A 
and B, we calculate intensities of the bremsstrahlung hard 
X-ray emission from the flaring loop. We use the simplified 


form of the bremsstrahlung cross-section (see e.g. Kontar 


et al. 2002): 


OO 

1(e) = const j J N(E,l) n ^dEdl. 


(26) 


The synthetic intensity maps are shown in Figures 15 
and [16] for four different moments. In order to compare 
them with observational data, they are smoothed by a 
Gaussian profile with the half-width of ~ 1.5 Mm (com¬ 
parable to the RHESSI spatial resolution). The main dif¬ 
ference between the low-density case (Model A) and high- 
density case (Model B) is that the latter shows noticeable 
extended emission from the loop, while the low-density case 
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Fig. 13. Normalised pitch-angle distributions of high-energy particles {£ > 5 keV) during fast reconnection stage (t = 
1120to)- Left panels (a and c) are for electrons; right panels (b and d) are for protons. Upper panels (a and b) correspond 
to the low-density case (Model A), while lower panels correspond to the high-density case (Model B). 


has noticeable emission only from the footpoints. This can 
be easily explained by comparing the electron mean-free- 
path in these two models: in the low-density model, the 
mean free path for the energies ~ 10 keV is substantially 
longer than the loop length (~ 20 Mm), while it is only 
~ 1 Mm in the high-density case. 


4. Summary 

In the present work, we use a combination of MHD and test- 
particle methods to study energy release and high-energy 
particle motion in a reconnecting twisted coronal loop. The 
approach is generally similar to that previously used in 2D 


I 


Gordovskyy et al. 


m 


.2011 ) and 3D models ( [Gordovskyy and 
rownmg} |ZUll[ 2012[ ): time-dependent electric and mag¬ 
netic fields ana density distributions are used as an input 
for guiding-centre test-particle calculations of proton and 
electron trajectories. The main benefits from the present 
model are that 


we use a realistic loop-like structure (similarly to Kliem 


et al. 2010), and the twist is created by slow footpomt 


rotation; 

— we take into account strong stratification in the lower at¬ 
mosphere and consider different densities in the corona; 

— we take into account the effect of Coulomb collisions on 
test-particle motion. 

Therefore, this study attempts to combine energy release 
during reconnection with particle acceleration and trans¬ 
port within a single model. Obviously, the main problem is 
that the feedback of energetic particle motion on the elec¬ 
tric and magnetic field evolution is ignored. However, this 
should not significantly affect evolution of the system in the 
cases studied here, since the acceleration efficiency is low: 
high energy particles carry less than ~ 10% and h ence , the 


plasma is almost “thermal” (see also Gordovskyy et al. 
20121 . 


Magnetic reconnection in the present models is simi¬ 
lar to that observed in cylindrical fluxtubes with analyt- 


ical magnetic field twist profiles (e.g. Hood et al. 2009 


Gordovskyy and Browning 20111. Thus, a strong, helically 


shaped current is formed after the kink occurs; the current 


stucture gradually becomes filamentary as the reconnection 
proceeds. The connectivity changes in two different ways: 
fir st, the magnetic twist is reduce d and second, as shown 
by Gordovskyy and Browning (2011), the twisted loop held 
lines reconnect with ambient held lines, resulting in an ef¬ 
fective increase in the loop cross-section. However, there 
are several features in the present model, which are not 
observed in models with cylindrical huxtubes. First, a cur¬ 
rent torus, almost perpendicular to the loop, is formed near 
the loop-top just after the kink. Second, although the mag¬ 
netic energy dissipates in three extended regions (near the 
footpoints and the loop-top), the held connectivity changes 
mainly due to magnetic reconnection just below the mag¬ 
netic torus, which is formed near the loop-top after the 
kink. Finally, dissipation of the azimutal held leads to a 
substantial shortening of the loop held lines. The density 
of the plasma in the corona does not seem to have a notice¬ 
able effect on reconnection, as expected, since the plasma 
beta is low (around 10 -1 ) in both cases. 

Particles can be accelerated almost anywhere within the 
reconnecting loop, although they are predominantly ener¬ 
gised close to the loop-top, where the ratio of electric held 
to Dreicer held is the highest. The energy spectra are com¬ 
binations of Maxwellian thermal distributions at low en¬ 
ergies (below keV) and are nearly power-low high-energy 
tails. The spectral indices for electrons in the 10 — 100 keV 
range are between 1.5 and 3.0, which means the spectra 
are slightly softer than in similar collisionless models. The 
electron spectra are harder in the case with dense coronal 
plasma (10 11 cnr 3 ) than in the case with normal coronal 
density (10 9 cm -3 ). The acceleration efficiency strongly de¬ 
pends on the ambient plasma density. This is not surpris¬ 
ing, since the particle acceleration generally depends on the 
ratio of electric field to the so-called critical electric field 
E C rit ~ mn/£ (i.e. the minimum electric field required to 
produce runaway particles; it is equal to th e Dreicer field 
in case of thermal electrons) (see Sect. 3.41. In both cases 
(low- and high-density corona), the electric field is higher 
than the corresponding E cr n value for electrons. However, 
in the high-density case the electric field in the system is 
comparable to the corresponding E cr n value for protons, 
and there is almost no proton acceleration. This results 
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Fig. 15. Synthetic hard X-ray emission at 10 keV from the simulated flaring loop in low-density atmosphere. 


in an interesting picture: in the low-density medium, the 
non-thermal particle energy is nearly equally split between 
protons and electrons, while in the high-density medium, 
protons remain nearly thermal amid a substantial number 
of non-thermal electrons. This can potentially explain why 
some flares effectively accelerate both electrons and ions, 
and other flares produce substantial number of electrons, 
while not showing any noticeable ion emission. Many phe¬ 
nomenological models suggest that this happens because, 
even though there are plenty of ions, their energy is insuffi¬ 


cient to produce nuclear de-excitation lines needed for de¬ 
tection. However, it might be the case that ions in electron- 
only flares do not get accelerated due to high plasma density 
in the acceleration region. 

High-energy particles are rather uniformly distributed 
along the flaring loop. There is a noticeable increase in the 
particle number towards footpoints, which can be easily 
explained by lower average parallel velocities due to mag¬ 
netic field convergence. There is a visible proton-electron 
asymmetry between footpoints at the early stages of ac- 
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Fig. 16. Synthetic hard X-ray emission at 10 keV from the simulated flaring loop in high-density atmosphere. 


celeration, because the initial configuration has a non-zero 
net current and hence, the electric field has a preferen¬ 
tial direction along the fluxtube before it gets filamentary. 
Although this proton-electron asymmetry would result in 
the development of an oppositely-directed large-scale elec¬ 
tric field along the reconnecting loop (the so-called return 
current electric field), this effect can still exist, as the non- 
thermal proton-electron current can be compensated by an 
ambient thermal electron current, provided the plasma is 
highly ionised and the number of high-energy particles is 


low (which is the case here) (see Zharkova and Gordovskyy 


2004 2005). In the high-density case, this effect is less no¬ 


ticeable - owing to pitch-angle scattering due to Coulomb 
collisions. 

In general, the structure of the synthetic hard X-ray 
sources is similar to that observed with RHESSI (e.g. Lin!, 


2011). There are high-intensity footpoint sources, which 
are flat (their height is only about 0.5 — 1.0 Mm), and 
weaker extended emission from the whole loop. The ratio 
of the footpoint and extended source intensities depends 
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on the plasma density in the loop. In the case of normal 
coronal density, the emission comes mostly from the foo- 
points, while the intensity of the extended source in the 
high-density case is comparable to the footpoint intensities. 
The size of the footpoint sources increases with time (due 
to expansion of the twisted loop as discussed in Sect. |3.5[ 
however, this effect is less prominent in hard X-ray due to 
the assumed temporal and spatial resolution of the instru¬ 
ment. 

Similar to Gordovskyy and Browning (2012), we assume 
that the appearance of a strong parallel electric field due 
to anomalous resistivity in fragmented current structures is 
the main mechanism that accelerates particles. Other mech¬ 
anisms, such as acceleration due to magnetic field contrac¬ 
tion and large scale MHD waves, are included within our 
model and can have an effect, although they are found to 
be negligible compared to the parallel electric field (for the 
assumed value of anomalous resistivity and the adopted 
length scale). 

Our work suggests that large-scale MHD turbulence de¬ 
velops during magnetic reconnection and results in frag¬ 


mentation of initially solid current (see also Browning et 


al. 2008; 


Hood et al. 2009). This phenomenon affects 
particle behaviour: proton and electron trajectories be¬ 
come more stochastic, as they are accelerated thr ough re- 
peated encounters wit h c urre nt frag ments (see Sect. 3.2 and 
Gordovskyy and Browning (2011 2012), see also Gar gill 


et al. (2012 [ ) for review). Another important factor, which 


can affect particle acceleration and transport is MHD tur- 


bulence at smaller spatial scales (e.g. 

Miller et al. 1996 

Larosa et al. 

1996; Cho and Lazarian 

2006 Cargill and 

Vlahos 

2009 

Bian and Kontar 

2010 

Bian et al.| 2010). 


connection in multiple small-scale curren t sheets (Lazarian 
and Vishniak 1999| Eyink et al. 2011). Hence, the pres¬ 


ence of small-scale turbulence, in principle, could drasti¬ 
cally change characteristics of the electric field in the model 
and lead to a different picture of particle acceleration. The 
implicatio ns of turbulence f or par ticle acceleration are sum¬ 
marized in Lazarian et al. (2012). 

Furthermore, turbulence can make substantial contribu¬ 
tion to particle scattering during their transport. Scattering 
due to Coulomb collisions considered in the present model 
is expected to be dominant in the dense chromosphere and 
transition region. At the same time, small-scale MHD tur¬ 
bulence can be important for particle scattering in the low- 
density corona, along with various types of self-induced 
kinetic waves. Thus, high-energy ions have been shown 


to generate kinetic Alfven waves (e.g. Voitenko 1998) 


which, in turn, result in noticeable heating in the corona 
(|Gordovskyy et al. 2005) (while collisional scattering re- 
sults only in effective heating of dense transition region and 
chromosphere). 

However, it is impossible to include these microscopic 
effects in our large-scale model. The macroscopic descrip¬ 
tion of these microscopic effects is required for them to be 
included in large-scale transport models, which necessitates 
more studies into this issue. 

The scenario developed in the present paper is a good 
candidate for interepreting small non-eruptive explosive 
events occuring in single loop configurations. A benefit of 
this model is reduced energy losses during particle trans¬ 
port due to more uniform (rather than localized at the 


loop-top) acceleration (see discussion in Gordovskyy et al. 


2012). The model, however, needs to be further improved; 


in particular, one needs to consider the possibility of con¬ 
tinuous helicity injection (or footpoint twisting). 
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